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We present a new and efficient method to obtain a Gamow shell-model basis and matrix ele- 
ments generated by realistic nucleon-nucleon interactions. We derive a self-consistent Hartree-Fock 
potential from the renormalized N 3 LO interaction model. The corresponding Gamow one-body 
eigenstates are generated in a plane wave basis in order to build a Gamow shell-model set of basis 
states for the closed shell nuclei 4 He and 16 0. We address also the problem of representing a realistic 
nucleon-nucleon interaction in a two-particle Berggren basis in the laboratory frame. To achieve this, 
an expansion of matrix elements of the residual nucleon-nucleon interaction in a finite set of har- 
monic oscillator wave functions is used. We show that all loosely bound and narrow resonant states 
converge surprisingly fast. Even broad resonances in these two-particle valence systems converge 
(N ' within a reasonable number of harmonic oscillator functions. Examples of 6 He and 18 O Gamow 

shell-model calculations using 4 He and 16 O as closed shell cores are presented. This procedure al- 
lows Gamow shell-model calculations to be performed with all realistic nucleon-nucleon interactions 
and with either momentum or position space representations for the Gamow basis. Perspectives for 
, nuclear structure calculations of dripline nuclei are outlined. 

r- ■ 

1 PACS numbers: 21.60.Cs, 21.10.-k, 24.10.Cn, 24.30. Gd 

(N : 

> ' I. INTRODUCTION 

in : 

A challenge in modern nuclear physics is the description of nuclei far from the valley of stability. These nuclei 
^vq ■ exhibit unusual features such as very low particle-emission thresholds, halo densities and unbound ground states. A 
(~) proper understanding of the mechanisms underlying the formation of such nuclei is presently a great challenge to 
nuclear theory, especially the case of two-neutron Borromean halos such as 6 He and 11 Li. The theoretical description 
of such exotic nuclei cannot be worked out within standard models because of the appearance of strong couplings to 
. the continuum. 

The extreme clusterization of Borromean nuclei into an ordinary core nucleus and a veil of halo nucleons has 
motivated few-body approaches such as the hyperspherical harmonic method and momentum space Faddeev equations 
to these nuclei pj. However, the few-body modeling of Borromean and halo nuclei is not completely satisfying as the 
treatment of core excitations and the anti-symmetrization between core and valence nucleons is therein approximate. 
An ab initio description of these nuclei, taking into account all relevant degrees of freedom, would alleviate the 
. ' defects of such cluster approaches. To achieve this, a reformulation of the shell model using a single-particle basis of 
bound, resonant and scattering states appears to be the most straightforward method. The Continuum shell model 
0- El [J IM an d the recently developed shell model embedded in the continuum (SMEC) 0, H, 0, offer such 
a possibility. In SMEC, two subspaces of bound/quasi-bound states and scattering states are introduced and their 
coupling taken into account following the techniques discussed in for example Ref. However, most calculations 

have been performed with only one-particle decay channels. While the theoretical formulation of SMEC with two- 
particle decay channels has been formulated (see Ref. [Toj with applications to two-proton radioactivity), exact 
three-body asymptotics have never been applied numerically. The very rapidly increasing complexity of SMEC with 
many-body decay channels is a hindrance to the study of cluster-emitting systems, such as in particular Borromean 
nuclei. 

The newly developed Gamow shell model [IlllllliQHHIllHHIll has proven to be a reliable tool 
in order to probe the structure of such nuclei. This model unifies structure and reaction properties of nuclei, and 
most importantly allows for an exact treatment of antisymmetry and has no limitation on the number of particles 
in the continuum. It is then particularly well suited for the study of Borromean nuclei. The starting point of the 
Gamow shell model is the B ergg ren completeness relation, where bound, resonant and scattering states are treated 
on an equal footing 0, |2^, |U 0> E3 ■ The completeness relation is built upon bound, resonant states and an 
integral over a continuum of scattering states with complex energy. This integral has to be discretized in order to be 
applied in numerical calculations. A complete many-body Berggren basis is then constructed with Slater determinants 
integrating bound, resonant and non-resonant discretized continuum orbitals. The Gamow shell model can be seen 
as a direct generalization of the standard shell model, where the standard harmonic oscillator set of states is replaced 
by a Gamow basis. 
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An important question concerns the choice of the potential to generate the one-body Gamow basis states. In Gamow 
shell-model calculations, the single-particle basis has normally been constructed from a Woods-Saxon or a Gaussian 
potential depicting 4 He or 16 cores, fitted to reproduce the single-particle states of 5 He and 17 0, respectively fillip . 
However, in a fully microscopic approach, the single-particle basis should be constructed from the free nucleon-nucleon 
interaction or more complicated three and/or many-body interactions. This can be done by summing various diagrams 
in many-body perturbation theory. At lowest order this approach is given by the Hartree-Fock approximation (see 
Ref. [Ilj where a Gamow-Hartree-Fock basis was derived and applied to schematic interactions.) 

In many-body perturbation theory, one cannot use the free nucleon-nucleon interaction, since it yields strongly 
repulsive and/or diverging matrix elements at short internucleonic distances. In order to remove these divergencies, 
renormalized nucleon-nucleon interactions have been constructed from the Brueckner G-matrix approach |26l l27l |28| . 
The G-matrix is a soft interaction, which is obtained by resumming in-medium particle-particle correlations. 

Recently, an alternative renormalization scheme which i nteg rates out the high momentum components of the 
nucleon-nucleon interactions has been proposed [2^, |3(], El l32l l33[. Using a similarity transformation of the two- 
nucleon Hamiltonian, a Hermitian soft-core effective nucleon-nucleon interaction is obtained in a model space defined 
by a cutoff A in the relative momentum between the nucleons. This renormalized interaction has become known as 
a low- momentum nucleon-nucleon interaction, labeled Vj ow -k- The interaction Vjow-k is an energy and nucleus inde- 
pendent effective interaction which reproduces nucleon-nucleon scattering data, but displays a sizeable dependence 
on A. 

In this work, the single-particle Gamow Hartree-Fock basis is constructed using a renormalized interaction of the 
Viow-k type, derived by similarity transformation techniques of the nucleon-nucleon interaction. Our renormalization 
scheme requires a plane wave basis formulation of the Schrodinger equation. Such a basis is a natural starting 
point since nucleon-nucleon interactions are usually derived explicitly in momentum space, as for example the N 3 LO 
interaction |34l l35l|. In order to perform Gamow Hartree-Fock calculations, the nucleon-nucleon interaction has to 
be defined by the coordinates of the laboratory system. The transformation of the interaction from the relative and 
center of mass frame to the laboratory frame is performed with the so-called vector brackets [13, HE Hll • These 
are the less known momentum space analogs of the Moshinsky transformation coefficients of the harmonic oscillator 
representation, generalizing the Talmi transformation to arbitrary bases. In the presence of unbound states such as 
in a Gamow basis, the single-particle potential has to be analytically continued in the complex fc-plane. In Ref. [H^| . 
it was shown how a single-particle Berggren basis can be obtained by the contour deformation method in a basis of 
spherical Bessel functions. 

For a microscopic approach to be fully consistent, the realistic nucleon-nucleon interaction should generate both a 
single-particle basis through Hartree-Fock calculations and an effective nucleon-nucleon interaction to be diagonalized 
in the Gamow shell model. We can obtain this by letting the renormalized nucleon-nucleon interaction to be ex- 
pressed in a two-particle Berggren basis. However, the difficulty in analytically continuing the vector transformation 
coefficients to the complex A-space, prevents such a derivation. In this work, an alternative approach to calculate 
realistic interactions in Gamow bases is proposed. The method is based on an expansion of the nucleon-nucleon 
interaction in a finite set of harmonic oscillator wave functions. Within this framework, the analytic continuation 
of the nuclear interaction is trivial, and matrix elements can therein be very efficiently calculated through the use 
of the standard Talmi transformation. As will be shown, this method provides well converged energies and wave 
functions in the Gamow shell-model calculations with a small number of harmonic oscillator states. In addition, as 
harmonic oscillator wave functions have a similar behavior in momentum and position space, both momentum space 
and coordinate space representations can be used for the Gamow basis. This method may also provide a solution to 
the problem of spurious center of mass motion in Gamow shell- model calculations. 

The outline of the paper is a follows. In Sec.[nj the derivation of a renormalized nucleon-nucleon interaction suitable 
for a perturbative many-body approach in the Gamow shell model is described. In Sec. lIIII the self-energy and Gamow 
Hartree-Fock single-particle basis of the Vi ow -k interaction are constructed and applied to the 4 He and 16 O closed-shell 
nuclei. Sec. II VI outlines the harmonic oscillator expansion method for the nucleon-nucleon interaction, and Sec. [V] 
illustrates applications in Gamow shell-model calculations for two selected valence systems, 6 He and 18 O. There we 
discuss also the convergence of narrow and broad resonances as functions of the number of harmonic oscillator wave 
functions used in the expansion. Sec. IVII points out the equivalence between the momentum and the position space 
formulations of the Gamow shell model when the harmonic oscillator expansion method is used. Finally, in Sec. IVIII 
we outline our conclusions and future perspectives. 

II. RENORMALIZED NUCLEON-NUCLEON INTERACTION 

In order to build the Gamow Hartree-Fock potential in fc-space and a Gamow shell-model Hamiltonian matrix, 
it is necessary to construct the self-energy E(fc a Z a j a , fc&) defined by the inclusion of various diagrams in many-body 
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perturbation theory, discussed in Sec. IIHI . Note here and in the following discussion the distinction between k a ,kb 
and k and l a and I. The notations k a or l a (latin letters) refer to the quantum numbers of a single-particle state a, 
whereas I or k without subscripts (or with greek letters as subscripts) refer to the coordinates of the relative motion. 

To compute many-body perturbation diagrams, the nucleon-nucleon interaction has to enter a perturbative treat- 
ment. Hence, the free nucleon-nucleon interaction, giving rise to diverging matrix elements, cannot be used directly 
and has to be renormalizcd. Since parts of our formalism is based on computing the self-energy in a momentum basis, 
it is convenient here to use a renormalization scheme based on a cutoff in momentum space as discussed by Bogner 
et al p9l ] and Fujii et al |30l l3l| . 

This approach is based on two steps, a diagonalization in momentum space for relative momenta k € [0, 00) of the 
two-body Schrodinger equation and a similarity transformation [30ll3]| to relative momenta k £ [0, A], A defining the 
relative momenta model space. Typical values of A are in the range of ~ 2 fm _1 . The nucleon-nucleon interaction is 
diagonal in the center of mass motion. One can therefore easily map the full diagonalization problem onto a smaller 
space via a similarity transformation and obtain thereby an effective interaction for a model space defined for low 
momenta. This interaction has been dubbed Viow-k in the literature, see for example Ref. j2£j- The effective low- 
momentum interaction Viow-k is constructed in such a way that it reproduces exactly the main characteristics of the 
nucleon-nucleon wave function in the full space. 

The interaction Vjow-k looks attractive at first glance, but may generate undesirable features in the Gamow shell 
model. Many-body calculations using a renormalized nucleon-nucleon interaction of the low-momentum type introduce 
a strong dependence on the cutoff A in momentum space. By integrating out high momentum modes of the nucleon- 
nucleon interaction, one excludes certain intermediate excitations in the many-body problem. While the two-body 
problem is exact with V[ow-k, the three-body problem will not be. In Ref. 33] Vi ow -k was accompanied with a A- 
dependent three-body force in order to reproduce the ground-state energies of 4 He and 3 H for each value of A. It is 
hoped that a three-body force is sufficient to eliminate the A dependence for heavier nuclei. However, if it turns out 
that one needs to go beyond three-body forces for nuclei with A > 4, many-body calculations starting with Vjow-k 
are futile. 

Alternatively, one could have defined a so-called G-matrix in momentum space as effective interaction p^.l27ll23 . l3^| . 
The latter introduces a dependence on the chosen starting energy and a reference to a given Fermi energy. This 
dependence can be eliminated by introducing for example higher-order terms in many-body perturbation theory |39| . 
We relegate such a discussion to future work. It must be stressed that the aim here is to demonstrate the feasibility of 
obtaining a single-particle basis for Gamow shell-model calculations using a realistic interaction. We adopt therefore 
a pragmatic approach and use Viow-k simply because it is easier to implement numerically in order to renormalize 
the nucleon-nucleon interaction. 

In the following we outline the procedure to obtain a Hcrmitian interaction Vi ow -k based on the similarity transfor- 
mation discussed in Refs. |30L l3ll I40L l41| . A unitary transformation can be parametrized in terms of the model space 
P and the excluded space Q via 

„_ f P(l + utu)- 1 / i P -Pwt(l + W«t)-1/2Q\ 

U ~ \ Qu(l + ^u)-y 9 P Q(1-|W)-1/2Q ), W 

where the wave operator oj is defined to satisfy the condition 

u = QojP, (2) 

the so-called decoupling condition Note that the unitary transformation is by no means unique. In fact, one can 
construct infinitely many different unitary transformations which decouple the P and the Q subspaces, as discussed 
by Kuo et al [43j. The above transformation depends only on the operator u) which mixes the P and Q subspaces 
and is in some sense "the minimal possible" unitary transformation. Following the method of Ref. |30j , one obtains 

U= (l+w-w^)(l+u)U^ + w^w)- 1/2 . (3) 

The above operator U leads to the effective interaction V using the definition 

V = U- 1 {T + V)U -T, (4) 

where T is the kinetic energy of the nucleons and V is the free nucleon-nucleon interaction. 

To express the renormalized interaction in momentum space, one starts with the Schrodinger equation for the 
relative momentum k, 

dk' k' 2 (k\T + V\k'){k'\^ a ) =E a (k\^j a ), (5) 
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where the plane wave states are eigenfunctions of the kinetic energy operator T in the relative system and form a 
complete set 

dkk 2 \k)(k\ = 1. (6) 

The momentum space Schrodinger equation is solved as a matrix equation by discretizing the integration interval by 
some suitable rule, (here, the Gauss-Legendre quadrature is used). The discretized Schrodinger equation reads 

^2w 7 k^{k s \T + V\k r ){k r \il> a ) =E a {k g \i> a ), (7) 

7 

where fc 7 are the integration points and w 7 the corresponding quadrature weights. Introducing \ks) — k$^/ws\ks) , 
Eq. J7J becomes 

Y,(ks\T + V\k 7 )(k^ a ) = E a (k s \^ a ), (8) 

7 

where 

N 

\h)(k s \ = 1, (fc# 7 ) = 8 5<1 . (9) 

5=1 

The matrix elements of the Hamiltonian are expressed in the plane wave basis : 

-ff<5, 7 = (ks\T + V\k 7 ) — —5s n + y/wsWjkskyViks, fc 7 ), (10) 

where m is the average of the proton and the neutron masses. The full space is now divided in a model space P and 
an orthogonal complement space Q. The model space P consists in the Np plane wave states lying below the cutoff 
A, and the Q-space consists of the remaining states, viz. 

P = {\k), \k\ < A} , Q = {\k), A < |fc| < oo} . (11) 

The model space is thus defined for all momenta k € [0, A] fm _1 . In order to obtain an effective interaction in the 
model space P, the decoupling condition in Eq. J5J has to be fulfilled. Once the transformation matrix ui in the plane 
wave basis (k\QuiP\k) is obtained, the low-momentum effective nucleon-nucleon interaction Vi ow -k reads 

(k\PV iow ^ k P\k') = ^^(fc|P(P + ^) 1 / 2 P|fc") (k"\P(T + V)P\k"') (k"'\P(P + o/VT 1/2 .P|F) 
k" k'" 
k 2 

Skk', (12) 

m 

see also Ref. for further details. The effective model space interaction in the original plane wave basis \ks) is then 
given by 

(fc<5|Mow-k|fc 7 ) = — , , — , (13) 

y/wsw-fksk-y 

where {\k s ), |fc 7 )} e P. 



III. DERIVATION OF SELF-CONSISTENT GAMOW HARTREE-FOCK BASIS 

The renormalized nucleon-nucleon interaction Mow-k is defined in terms of various quantum numbers as follows 

(klKL(J)ST z \ T/ low _ k \k'l'KL(J)ST z ) , (14) 

where the variables k, k' and I, I' denote respectively relative and angular momenta, while K and L are the quantum 
numbers of the center of mass motion. J , S and T z represent the total angular momentum in the relative and center 
of mass system, spin and isospin projections, respectively. 
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The A— body Hamiltonian H is defined as 

A A 

ff=— $>? + ^ ow _ k (i,j)- (!5) 
2m ' 

l=l Kj 

The spurious center of mass energy is removed by writing the internal kinetic energy as 

v ' 8=1 i<j 

The introduction of an additional two-body term yields a modified two-body interaction 

A , k k \ 

Hi = Mow-k + K, n . = ( tW-k(*,i) - ) , (17) 



resulting in a total Hamiltonian given by 



A 



v ' i=i 

In all calculations reported in this manuscript, the modified two-body Hamiltonian defined in Eq. (|17|l is employed. 

Starting from the renormalized momentum-space version of the nucleon-nucleon interaction Vj ow -k, with matrix 
elements in the relative and center of mass system system, one can obtain the corresponding matrix elements in the 
laboratory system through appropriate transformation coefficients |36l 13^. l38| . This transformation proceeds through 
the definition of a two-particle state in the laboratory system. With these coefficients, the expression for a two-body 
wave function in momentum space using the laboratory coordinates can be written as 



\(kJajat Za )(k b l b J b t Zb )JT z ) = J k2<lk I R2dK 

ILXSJ 




(19) 



x(-iy +J - L - s Fj\*jj b s 



L I A 
S J J 



x (klKL\k a l a k b l b X) \klKL(J)SJT z ) , 

where (klK L\k a l a k b l b \) is the transformation coefficient (vector bracket) from the relative and center of mass system 
to the laboratory system defined in Refs. [13, |3^. The factor F is defined as F = (1 — (— l) l+s+Tz ) / \/2 if we have 
identical particles only (T z = ±1) and F — 1 for different particl es (p rotons and neutrons here, T z = 0). 
The transformation coefficient {klKL\k a l a k b l b ) is given by [36l l37l l3Sl | 

4?r 2 

(klKL\k a l a k b l b X) = — 5(w)0(l - x 2 )A(x), (20) 



with 



w = k 2 + \K 2 -±(k 2 a + k 2 b ), (21) 



x = {kl - k 2 - -K 2 )/kK, 



(22) 



and 



A(x) = 



2A + 1 ^ 



^<Y l (k) x Y L (K)] X * x \Y l >(k a ) x Y l »(k b )ti. 



(23) 
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(a) (b) (c) 

FIG. 1: Diagrams through second order of the interaction Hi. Diagram (a) is the first order Hartree-Fock term while diagrams 
(b) and (c) are the 2p — lh and 2h — lp corrections, respectively. The labels a and c represent the outgoing and incoming states, 
respectively. The intermediate particle states are labeled as p,p\,p2 with upward arrows and the intermediate hole states as 
h, hi, fi2 with downward arrows. The dotted lines represent the interaction Hi. 

The functions Y l are the standard spherical harmonics, x is the cosine of the angle between k and K so that the step 
function takes input values from to 1. In our codes, the coordinate system of Kuo et al. (3^ was chosen. 

To compute the Hartree-Fock (HF) diagram (see Fig.^(a)) we need matrix elements in a mixed representation of 
bound and scattering states such as 

({kJajat Za ){n b l b j b t Zb )JT z \ Hi \{k c l c j c t Za )(n d l d j d t Zd )JT z ) , (24) 

where the labels a and c represent scattering states and b and d represent bound states. All matrix elements discussed 
here are assumed to be antisymmetrized (AS). Note that the two-body center of mass correction term V c . m . in Eq. i|17|) 
is calculated directly in the laboratory coordinates using the Wigner-Eckart theorem, while the renormalized nucleon- 
nucleon interaction is given in laboratory coordinates using the transformation given in Eq. I|19|l . The calculation of 
these matrix elements requires the knowledge of two-body states in a mixed representation with for example harmonic 
oscillator wave functions R na i a representing the bound states and plane waves for the resonant or continuum states 

\{n a l a j a t Za )(k b l b j b t Zb )JT z ) = J k 2 a dk a R na i a {k a )\{k a l a j a t Za ){k b l b j b t Zb )JT z ) , (25) 

where k a l a j a and n a l a j a are respectively plane wave and harmonic oscillator wave functions. The two-body state is 
represented by the quantum numbers of the total angular momentum J and isospin projection T z . 
With these matrix elements, the expression for the HF diagram shown in Fig. ^a) can be derived 



VHF{k a k c l a j a t Za ) = — 51 J 2 ((kJajat Za )(n h l h j h t Zh )JT z \ Hi \(k c l a j a t Za )(n h l h j h t Zh )JT z ) AS , (26) 

Ja J n h l h j h tz h 

where x = y/2x + 1 and nhlhjht Zh are the quantum numbers of the nucleon hole states. The variables l a , ja, t Za are the 
orbital angular momentum, total angular momentum and isospin projection (t Za = ±1/2) of the incoming/outgoing 
nucleon, and k a (fc c ) the outgoing (incoming) particle momenta. The mixed matrix elements of the two-body center 
of mass correction term needed in the HF calculation reads 

{{kalajat Za ){n h l h j h t Zh )JT z \ - ^ ^ 2 \(k c l c j c t Zc )(n h l h j h t Zh )JT z ) AS = 
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( 1 + (-l) l «+ lh+1 ] 

X | |*a*c-R nfc I fc (*a)-R« h i h (*!c)«i. J - *I a J.<5t..*. h *t..t. h - ( 2 ?) 

In summary, if we limit ourselves to the computation of the HF contribution, the expression for the self-energy 
reads 

^(jJakakc) =VHF(jJak a k c )- (28) 

In this work, only the HF contribution is considered, while in Ref. [2(| the authors also studied contributions from 
2p — lh and 2h — lp intermediate states. They yield an imaginary term which can be related to its real part via a 
dispersion relation. 

To calculate the contributions from the 2p—lh diagrams like the example displayed in Fig. [2(b) (or similarly the 
2h — lp diagram of Fig.^(c)) we evaluate the imaginary part first. The real part is obtained through the dispersion 
relations to be defined below. The analytical expression for the imaginary contribution of the 2p—lh diagram, which 
gives rise to an explicit energy dependence of the self-energy, is 



m P -ih(jJak a k c t Zci cj) = Yl Z Z / k2dk J K 2 dKJT 

3a n h l h j h t Zh J ILSJ 

x {(kJajat^irihlhjht^JT^Hi \klKL(J)SJT z ) 
x (MKL(J)SJT Z \ Hi \(k c l a j a t Za )(n h l hJh t Zh )JT z ) 

x *\" +£h -m- N -M- N ) < (29) 

where u> is the energy of the incoming nucleon in a state a. The quantities klKL{^f)SJT z are the quantum numbers of 
the intermediate two-particle state. To compute the two-particle-one-hole diagram given by the second-order diagram 
of Fig.n(b), the following matrix elements are needed 

({kJ a j a t z J{n b l b j b t Zb )JT z \Hi \klKL{J)ST z ) . (30) 

The contributions to the real part of the self-energy from Eq. I|29l) can be obtained through the following dispersion 
relation 

, , , , v P f°° yV2 P -ih(jJ a k a k c t Za uj') , 

V2p-lh{jJakakct Za LU) = — dlV , (31) 



where P takes the principal value of the integral. Since W2 P -ih is different from zero only for positive values of ui' and 
its diagonal matrix elements are negative, this dispersion relation implies that the diagonal elements of V2 P -ih will 
be attractive for negative values of to. This attraction should increase for small positive energies. It will eventually 
decrease and become repulsive only for large positive values of the energy of the interacting nucleon. Similar expressions 
can also be derived for second-order diagrams with 2h — lp intermediate states 26] . Inclusion of these terms will be 
presented in a future work (see also the discussion in Sec. IVIlj l. 

The equations above for the nucleon self-energy are only valid along the real-energy axis. However, ~E(j a l a k a k c ) has 
to be analytically continued from the real fc-axis to the complex fc-plane in order to obtain a genuine Gamow shell- 
model single-particle basis. The derivation of a Berggren basis in momentum space with the Contour Deformation 
Method was described in Ref. 19] . The analytical continuation of T,(j a l a k a k c ) to the complex fc-plane leads to 
theoretical and practical difficulties due to the appearance of Dirac and Heaviside distributions in vector brackets (see 
Eq. H20fl ). However, once the self-consistent self-energy has been obtained along the real fc-axis, the HF potential can 
be simply continued to the complex plane via two sets of Fourier-Bessel transformations. To obtain a self-consistent 
HF potential in the complex fc-plane, the following scheme is employed : 

• The first step is to calculate Vhf self-consistently on a grid on the real momentum axis using the interaction 
Hi of Eqs. l(T7|) and lt^|l . 

• When a self-consistent solution has been obtained, Vhf is calculated in position space via a double Fourier-Bessel 
transform. 

c\ poo roo 

VHF(jlrr') = - / dkk 2 / dkk 2 ji{kr)ji{k'r')V H F{jlkk'). (32) 



o 
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FIG. 2: Convergence of the Op-j/2 and the Opi/2 energies in 16 with respect to iteration number in the self-consistent HF 
calculation. Here we used a model space A = l.Ofm" 1 in the construction of Vi ow -k. 



• Having obtained Vhf in the r-plane we may go back to the complex fc-plane using one more Fourier-Bessel 
transformation. On an inversion symmetric contour L + in the complex fc-plane, the HF potential becomes 

/>oo />oo 

V HF (jlkk') = / drr 2 / dr'r' 2 j i (fcr)j i (fcV')V ffF 07rr'), (33) 
Jo Jo 

where fc and fc' belong to the contour L + and are therefore complex. The analytically continued single-particle 
Schrodinger equation on a general inversion symmetric contour then takes the form 

k 2 ipnlj{k)+ dk'k' 2 VH F (jlkk')^ n i(k') = E nl ip nl (k), (34) 

m eS J L + 

with m e g = 2to(1 — 1/A)^ 1 . Here both fc and fc' are defined on an inversion symmetric contour L + in the lower 
half complex fc-plane, resulting in a closed integral equation. In order to solve this equation, the integral has to be 
discretized, and we finally end up with a complex symmetric matrix diagonalization problem, in analogy with Eq. J7J|. 
This procedure results in a self-consistent Gamow Hartree-Fock basis, which is complete within the discretization 
space, and includes bound, resonant and a finite set of non-resonant continuum states (see Ref. [l^ for details about 
one-body Berggren completeness relations in momentum space). 

Based on this approach, the single-particle Gamow Hartree-Fock states of the closed-shell 4 He and 16 O are calculated 
with a low-momentum nucleon-nucleon interaction constructed from the realistic N 3 LO nucleon-nucleon interaction 
[34l l35|. For 16 we used the models spaces defined by A — 1.9, 2.0 and 2.1 fin -1 , and for 4 He the models spaces 
defined by A = 1.8, 1.9 and 2.0 far 1 

Table [I] presents the neutron single-particle energies obtained with a 16 O core. We note that holes states are in 
general overbound for all chosen cutoffs A and that the spin-orbit splitting between the p 3 / 2 and p 1 / 2 states is too 
large. The particle states are in better agreement with data and for all model spaces the d 3 / 2 state is the only one 
which comes out as a resonance, in agreement with experiment. The spin-orbit spacing the d§/ 2 and d 3 / 2 states is 
also fairly well reproduced. Fig. [21 shows the convergence of the 0p 3 / 2 and the Opi/ 2 spin-orbit partners in 16 O with 
respect to the number of iterations in the self-consistent HF calculation. 

Table HH presents the calculated self-consistent neutron single-particle energies with respect to a 4 He closed shell 
core. The calculated values for the p 3 / 2 and the P1/2 energies are not so far from the experimental values. For a model 
space A = 1.8 fm _1 , the obtained width of the p 3 i 2 resonance coincides with the experimental width 0.648 MeV, 
while the calculated width of the pi/ 2 resonance (~ 7.4 MeV) is larger than the experimental value 5.57 MeV. The 
spin-orbit splitting between the P3/ 2 and the Pi/ 2 levels is fairly well reproduced, the experimental value is 1.27 MeV 
|45| while our values vary from 1.36 to 1.68 MeV. Noting that our calculations are done at the Hartree-Fock level, 
there is clearly room for improvements. However, although our results for hole states are overbound, we obtain a 
qualitatively correct spectrum. Higher-order corrections should improve the agreement with experiment, since 2p — lh 
and 2h — Ip correlations provide additional binding and improved spin-orbit splittings for particle states. 

Many-body calculations using a renormalized nucleon-nucleon interaction of the low-momentum type will unfortu- 
nately introduce a model space dependence in momentum space (see Tables III and IhTi . The model space dependence 
can only be eliminated by introducing the corresponding three- and many-body forces which the low-momentum two- 
body interaction induces. In Gamow shell-model calculations, the inclusion of three-body forces is not feasible for the 
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TABLE I: HF calculation of single-particle energies in 16 using the low-momentum N 3 LO nucleon-nucleon interaction for 
three different model spaces. The single-particle energies E are given in MeV for both real and imaginary parts. Experimental 
data are from Ref. (44). 





A = 1.9 fm _1 


A = 2 mi" 1 


A = 2.1 fin -1 


Expt. 




lj 


Re[E] 


Im[E] 


Re[E] 


Im[E] 


Re[E] 


Im[E] 


Re[E] 


Im[E] 


Sl/2 


-73.977 


0.000 


-68.496 


0.000 


-63.078 


0.000 


-44.000 


0.000 


P3/2 


-37.082 


0.000 


-33.824 


0.000 


-30.650 


0.000 


-21.840 


0.000 


Pl/2 


-23.981 


0.000 


-21.749 


0.000 


-19.680 


0.000 


-15.664 


0.000 


^5/2 


-5.060 


0.000 


-3.810 


0.000 


-2.541 


0.000 


-4.143 


0.000 


Sl/2 


-3.531 


0.000 


-2.556 


0.000 


-1.781 


0.000 


-3.273 


0.000 


^3/2 


5.189 


-1.669 


5.353 


-1.928 


5.419 


-2.155 


0.937 


-0.048 



TABLE II: Same as in Tab. Q for 4 He. Experimental data from Ref. @. 

A = l.Sfm' 1 A = l.gfm" 1 A = 2.0fm~ 1 Expt. 

lj Re[E] Im[E] Rc[E] Im[E] Rc[E] Im[E] Re[E] Im[E] 
si/a -25.731 0.000 -24.541 0.000 -23.079 0.000 -20.578 0.000 
p 3/2 0.819 -0.325 1.041 -0.479 1.287 -0.667 0.890 -0.324 
Pi/2 2.497 -3.697 2.514 -3.777 2.648 -4.029 2.160 -2.785 

moment. If one wants to minimize the effect from many-body forces, a better approach might be to use a G-matrix 
for nuclear matter as in Refs. [26ll39| . However, higher-order correlations such as 2p— lh or 2h— lp contributions are 
necessary in order to minimize the dependence on the starting energy and the chosen Fermi energy |39| . Whether the 
model space dependence of Vjow-k can be softened by the inclusion of higher-order correlations such as 2p — lh and 
2h — lp will be investigated in a future work. Furthermore, although the G-matrix carries an explicit starting energy 
dependence, this dependence is needed when one wants to compute for example spectral functions. 

IV. MATRIX ELEMENTS OF REALISTIC NUCLEON-NUCLEON INTERACTION WITH A GAMOW 

HARTREE-FOCK BASIS 

As discussed in the previous section, starting from matrix elements of the nucleon-nucleon interaction in the relative 
and center of mass system system, one can obtain the correspondin g m atrix elements in the laboratory system through 
appropriate transformation coefficients, see for example Refs. |36l l37l l38| | and the discussion in the previous section. 
This transformation proceeds through the definition of a two-particle state in the laboratory system using the vector 
bracket transformation. However, the latter is very complicated to handle in practical calculations beyond the Hartree- 
Fock level. One has also to face the problem of two-particle intermediate states not orthogonal to the incoming and 
outgoing states (see Eq. (I29|) )■ Another method is then needed in order to efficiently calculate effective nucleon-nucleon 
matrix elements in the complex fc-plane and laboratory frame with a Gamow shell-model basis. 

The renormalized nucleon-nucleon interaction in an arbitrary two-particle basis in the laboratory frame is given by 

(a6|V[ow-k|cc;) = ((n a l a j a t Za )(n b l b j b t Zb )JT z \ Vi ow _ k \{n c l c j c t Zr ){n d l d j d t Zd )JT z ) . (35) 

The two-body state \ab) is implicitly coupled to good angular momentum J. The labels n a ...d number all bound, 
resonant and discretized scattering states with orbital and angular momenta (l a ...d, ja...d)- 

In order to efficiently calculate the matrix elements of Eq. (|35H . we introduce a two-particle harmonic oscillator 
basis completeness relation 

J2\aP)(ap\ = I, (36) 

a</3 

where the sum is not restricted in the neutron-proton case. We introduce the greek single particle labels a, (3 for the 
single-particle harmonic oscillator states in order to distinguish them from the latin single-particle labels a, b referring 
to Gamow states. The interaction can then be expressed in the complete basis of Eq. (|36[1 

v °" = EE |a/3)(a/3|Mow-k|7<5)(7<5|, (37) 

a<(3 7<i 
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where the sums over two-particle harmonic oscillator states are infinite. The expansion coefficients 

(a/3|Mow-k|7<5) = {{n a l a j a t Za )(nplpjpt Zl3 )JT z \ VLw-k \(n~ / l- ( j~ ( t Zi )(nskj8t zs )JT z ) , (38) 

represent the nucleon-nucleon interaction in an antisymmetrized two-particle harmonic oscillator basis, and may easily 
be calculated using the well known Moshinsky transformation coefficients, see for example Ref. for expressions. 

Matrix elements of Eq. I|35[) are calculated numerically in an arbitrary two-particle Gamow basis by truncating the 
completeness expansion of Eq. (|37|) up to N harmonic oscillator two-body states 

N N 

(ab\V osc \cd) « Yl £<a6|a/3}<a/3|tW-k|7*><7*|cd>. (39) 

a<f3f<S 

The two-particle overlap integrals (ab\af3) read 

{ab \ap) = mmzt^r^im (40 ) 

for identical particles (proton-proton or neutron-neutron states) and 

(ab\a(3) = {a\a){b\p) (41) 
for the proton-neutron case. The one-body overlaps (a\a) are given by, 

(a\a)= f dTT 2 <p a (T)R a (T)5 u J 3aJ J tata , (42) 



where ip a is a Gamow single-particle state and R a is a harmonic oscillator wave function. The radial integral is either 
evaluated in momentum or position space, indicated by the integration variable r. The important point to notice 
is that the only numerical calculations involving Gamow states are the overlaps of Eq. 1|42|1 . Hence, this expansion 
provides a simple analytical continuation of the nuclear interaction in the complex fc-plane. More precisely, the 
expansion coefficients of Eq. (|38|) can always be calculated with real harmonic oscillator wave functions, and in the 
case of Gamow functions spanning a complex contour L + (as in the momentum space representation), it is only the 
overlap integrals of Eq. (|42[1 which are analytically continued in the complex plane. These one-body overlap integrals 
converge in all regions of the complex plane which are of physical importance, due to the gaussian fall off of the 
harmonic oscillator wave functions. 

The convergence with the number of harmonic oscillator states N of the nuclear interaction expansion of Eq. (|38() can 
however not be checked by considering the matrix elements of Eq. (|39|l when N — > +oo. Indeed, they generally diverge 
when N — > +oo, due to the long-range character of the nuclear interaction in laboratory coordinates. This reflects 
the fact that the convergence of the Gamow shell-model Hamiltonian with N is weak, because the representation of 
a Hamiltonian such as Vjow-k in terms of a continuous Gamow basis is a distribution. Actually, only eigenvalues and 
eigenfunctions of the Gamow shell-model Hamiltonian converge with N, which will be shown in particular cases in 
Sec. E| 



V. GAMOW SHELL-MODEL CALCULATIONS OF 6 HE AND ls O. 



In Sec. IIIII we constructed a self-consistent single-particle Gamow Hartree-Fock basis starting from a realistic 
renormalized nucleon-nucleon interaction. Ultimately this basis should be used to construct the effective interaction 
in a given model space and diagonalize the shell-model Hamiltonian. Using the harmonic oscillator expansion method 
of the nucleon-nucleon interaction, we have a practical way of constructing the effective interaction to be incorporated 
in Gamow shell-model calculations. It is our aim to investigate whether a finite truncation of the harmonic oscillator 
expansion given in Eq. (|39|l may yield converged energies and wave functions in Gamow shell-model calculations. As 
a first application we consider two nucleons moving outside a closed core. The Gamow shell-model Hamiltonian used 
reads 

H(l, 2) = & HF (1) + ^hf(2) + "Ml, 2), (43) 

where /ihf is the self-consistently derived single-particle Hartree-Fock potential and V c s (1, 2) is the effective interaction 
acting between valence particles. In this work, all calculations are implemented up to first order in many-body 
perturbation theory, where the effective interaction is defined as 

U cff (l,2) = Mow-k(l,2)-^^. (44) 
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Im[k] 




FIG. 3: Contour L + in the complex fc-plane used in construction of the single-particle Berggren basis. The contour in specified 
by the points A,B and C discussed in the text. 



It should be noted that we utilize the harmonic oscillator expansion method only for the nucleon-nucleon interaction 
part, i.e. T4 ow -k(l, 2) rs V osc , while the center of mass correction term is treated exactly in the laboratory frame. The 
matrix elements of the center of mass term are given by 

({nJajatz a ){nblbjbt Zb )JT z \'k 1 ■ k 2 \{n c l c j c t Zc )(n d l d j d t Zd )JT z ) = 

-h 2 (-iy°+1«+ J I Ja Jb J \{nJ a j a \\Vi\\nJ c j c ){n b l b j b \\V 2 \\n d l d j d ), (45) 
[3d Jc 1 J 

where the expression for the reduced matrix elements can be found in Ref. 46]. In our Gamow shell-model study of 
the two-particle valence systems with realistic interactions, 6 He and 18 0, a two-particle model space built from the 
s i/2, P3/2, Pi/2, d 5 /2 and d 3 / 2 single-particle states is used. Each combination of the quantum numbers Ij consists of 
25 single-particle orbitals, totalling 115 orbitals. The same integration contour L + in the complex fc-plane is used for 
all partial waves (see Fig.|3J. In our 6 Hc calculations the contour is defined with A — 0.28 — 0.12i fm -1 , B = 0.5 fm -1 
and C = 4 fm' 1 , and in our 18 calculations it is defined with A = 0.52 - 0.12i fm _1 ,i? = 0.64 fm -1 and C = 4 
fm _1 . Using Gauss-Legendre quadrature, the discretization of L + has been carried out with 5 points in the interval 
(0, .A), 7 points in the interval (A, B) and 13 points in the interval (B,C). Convergence of the two-particle states 
with respect to the number of integration points has been checked. Our discretization of L + yields a precision of the 
energy calculation better than 0.1 keV for all states considered. In the following discussion a model space defined by 
A = 1.9 fm -1 is employed. The oscillator length is fixed at b = 2 fm. 



A. He results 



Table Hill gives the convergence of the Of ground state and the 2 J excited state of 6 He as functions of increasing 
number of nodes in the harmonic oscillator expansion of the interaction. A remarkable observation is that the Of 
ground and the 2f excited states of 6 He converge rather fast with respect to the number of harmonic oscillator 
functions, since n max = 10 is sufficient to reach convergence. Our calculations are comparable with the experimental 
values of —0.98 MeV for the + ground state and 1.8 — 0.06i MeV for the 2 + excited state in 6 He. Especially the 
2 + excited state is well reproduced. A splitting of ~ 1.5 MeV is obtained between the + ground state and the 2 + 
excited state, to be compared with the experimental value of 1.8 MeV. The binding energy of the + ground-state 
may be improved by going beyond first order in many-body perturbation theory, as shown in for example Ref. |28| . 
The 2p — lh and Ip — 2h diagrams may yield extra binding and improve spin-orbit splitting for the HF single- 
particle states of 4 He. In addition, it is well-known that the two-body core-polarization contributions improve the 
spectroscopy of systems with two and more valence nucleons |28l . These findings agree also with Brueckner-Hartree- 
Fock calculations for oxygen isotopes, see for examples Ref. 47] . At first order in perturbation theory, the spectrum is 
very much compressed. The agreement with experiment is partly improved with the introduction of core-polarization 
contributions. 

It can be concluded that the energies considered here for 6 He, converge with respect to the number of harmonic 
oscillator functions in the expansion of the nucleon-nucleon interaction. However, this does not imply that other 
observables converge with the same speed, especially those which are sensitive to the tail of the wave function. 
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TABLE III: Convergence of the 0i + and the 2i + energies in 6 He as functions of the number of harmonic oscillator nodes in 
the expansion of the realistic low-momentum N 3 LO nucleon-nucleon interaction. The model space is given by A = 1.9 fm -1 . 
The harmonic oscillator length is chosen at b — 2 fm. Energies are given in units of MeV. 
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FIG. 4: Plot of the radial density of the + ground state of He for three different harmonic oscillator expansions with 
respectively n max = 4, 10 and 16. 



In order to investigate such a dependence, the single-particle radial density operator is considered for the Of ground 
state in 6 He. The single-particle radial density operator is given by 

AT 

P = J E\n)(r i \, (46) 

i 

where TV is the total number of valence particles (N = 2 for 6 He in a 4 He core). This operator measures the 
probability for either particle 1 or 2 is to be found at the position r^. Fig. 0] shows plots of the diagonal part of p, that 
is p(r = ?*i = r2) for the + ground state of 6 He, with n max = 4, 10 and 16 in the harmonic oscillator expansion of the 
nucleon-nucleon interaction. In Fig. 21 there is no observable difference between the n max = 4, 10 and 16 results. In 
Fig. [3 we examine the tail of the wave function. The densities obtained for ?i max = 10 and 16 are indistinguishable, 
while the n max = 4 displays a very small deviation from the converged results. 



B. 18 O results 



Tab. IIVI gives the convergence of the 0i + ,02 + ,4i + and 42 + state energies and Tab. |V] the convergence of the 
2i + ,22 + ,23 + and 24 + of 18 0, as the number of nodes in the harmonic oscillator expansion increases. All of the 
states converge with a reasonable low number of harmonic oscillator functions n max ~ 10. Our calculation of the 
Of ground-state energy comes at — 12.23MeV which is very close to the experimental value of — 12.18MeV. The 
calculated splitting between the Of ground-state and the Of excited state is ~ 3.73MeV which is also very close to the 
experimental value 3.63MeV. We are also able to predict that the first excited resonant state is the 42 + state coming 
at energy —1.44 — 0.74iMeV, which is in agreement with experiment. However, we are not able to correctly describe 
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the splitting between the 0i + ,2i + and the 4i + states. In our calculations the 0i + and 2i + are almost degenerate. 
The discrepancy with experimental data is expected to be reduced by going beyond first order in perturbation theory, 
including contributions such as the core-polarization diagrams to the effective interaction. In Ref. |2f| , it was shown 
how the splitting of the 0i + , 2i + and the 4i + states of 18 is indeed improved by going to higher order in perturbation 
theory. In order to improve our Gamow shell-model calculations for 18 starting from realistic interactions, we must 
include higher order diagrams in the perturbation series for the nucleon self-energy and the effective interaction. This 
is a topic which will be followed up in the future. 

TABLE IV: Convergence of the 0i + ,02 + ,4i + and the 42 + energies in 18 as functions of the number of harmonic oscillator 
nodes in the harmonic oscillator expansion. A model space defined by A = 1.9fm was used. The harmonic oscillator length 
is fixed at b = 2 fm. Energies are in units of MeV. 
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VI. EQUIVALENCE BETWEEN POSITION AND MOMENTUM SPACE REPRESENTATIONS 

While the momentum representation of one-body Gamow states has been used in 0,113] and here, the Gamow shell- 
model was first introduced employing a position space representation, see for example Refs. pL4L Il6j|. A momentum 
space representation has however normally been preferred in constructions of effective interactions based on realistic 
nucleon-nucleon interaction models. There are several reasons for this. Realistic interactions are usually derived in 
/c-space, so that plane wave expansions are the most natural bases. Moreover, in connection with Gamow shell-model 
calculations, the momentum space representation was meant to lead to a faster convergence with the number of 
discretized scattering states, as Gamow wave functions in momentum space are usually more localized compared to 
those in the coordinate representation. Moreover, the momentum space representation of the one-body Schrodinger 
equation is an integral equation in the general case, contrary to its integro-differential form in r-space, known to 
be much more difficult to solve. The necessity of imposing the asymptotics of the Gamow state in r-space was also 
thought to give rise to a slower convergence with the number of discretized scattering states. Finally, the use of 
complex scaling to calculate two-body matrix elements in the position representation leads to extremely slow 
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TABLE V: Convergence of the 2i + ,22 + ,23 + and the 2i + energies in O as functions of the number of harmonic oscillator 
nodes in the harmonic oscillator expansion. A model space defined by A = 1.9£m _ was used. The harmonic oscillator length 
is fixed at b — 2 fm. Energies are in units of MeV. 
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calculations and cannot even regularize a large class of infinite matrix elements occurring in long-range interactions. 
It was this last nuisance which motivated the use of surface-peaked interactions in position space calculations, namely 
the Surface Delta Interaction (SDI) [l4lll5j and the Surface Gaussian Interaction (SGI) [TlLll^ . One will see however 
that both representations are in fact equivalent theoretically but also numerically, where the computational cost to 
obtain a given precision is comparable in both cases. The main point is that the possibility to use the harmonic 
oscillator expansion method of Sec. II VI removes all the problems previously encountered with both representations. 

In order to obtain a self-consistent Gamow Hartree-Fock basis in position space, one has to solve the one-body 
integro-differential Schrodinger equation 

<y(r) + / dr' V HF (jlrr')v nlj {r') = E nl(fnlj (r), (47) 

where Vjjf (j'Zrr') is the self-consistent HF potential given in Eq. Q32[). The wave function \<fnij) has to exhibit a 
pure outgoing wave function behavior for bound and resonant states for r — > +oo, whereas it has both incoming and 
outgoing components if it is a scattering state. In the last case, n must be understood as representing its wave number 
fc. Even though it is an integro-differential equation, it can be solved with standard methods with the use of locally 
equivalent potentials 0> so that its integration has the same complexity as differential equations occurring with 
purely local potentials. This generates a self-consistent procedure, as the locally equivalent potential depends on the 
state that it generates. This is of no importance in practical situations, as the integrated potential will be non-local 
in r-space only if it is a HF potential generated by a finite-range interaction. As HF pot entials have to be solved 
self-consistently, no numerical overhead can occur. The slow convergence noted in Ref. [14J was due only to the use 
of the trapezoidal rule to discretized the non-resonant continuum. With the use of the Gauss-Legendre integration, 
as performed in fc-space calculations from the beginning, results have improved dramatically in r-space calculations, 
reaching fc-space calculations quality |Hflj . 

In fact, the fundamental difference between the r and the fc representations for Gamow shell- model applications 
lies in their different discretization schemes. In fc-space, it is the Bessel completeness relation of Eq. (JBJ which is 
discretized. The \ij>nij) states are then obtained by diagonalization of VHF(jlkk') in the discretized Fourier-Bessel 
basis space. They will be denoted as \ip n i) Dk - In r-space, it is the completeness relation spanned by the \<p n i) states 
themselves which is discretized. Indeed, one has : 




AT 

^2w n \<p n i){tp„i\ ~ 1, (49) 
i=l 

where n 6 (6, d) means that one sums over all bound (b) and decaying (d) states above the contour L+. The first 
completeness relation, exact, becomes the second discretized completeness relation, approximate, where w n is 1 for 
bound and resonant states, and the Gauss-Legendre weight for scattering states. As a consequence, the \tp n i) states 
of r-space, that we label \<f n i) Dr are exact up to numerical precision, since they come from a direct integration of the 
Schrodinger equation. Approximations arise only from their discrete and finite number in Eq. (|49|l . The corresponding 
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\4'ni) Dk states have to be approximate, as they are generated by a finite number of Bessel basis states, meaning that 
there are not enough frequencies to expand them exactly. As a consequence, one has \tp n i) Dk — > \tp n i) Dr only at the 
continuum limit, that is for N — > +00. 

This difference may appear for observables that depend on large values of r or k, such as particle densities at large r 
or momentum densities at large values of k. But this is of no importance as discretization effects become preponderant 
in these regions, and may thereby most likely lead to numerically unstable results. 

Consequently, both representations can be used in shell-model problems without any loss of precision. The re- 
maining question of the method to handle two-body matrix elements in purely numerical bases has been answered in 
Sec. II VI with the use of harmonic oscillator expansions. One has seen in the latter section that the only numerical 
calculations involving Gamow states are the overlaps between Gamow and harmonic oscillator states of Eq. H42fl . 
which are obviously fast numerically. The Gaussian decrease of harmonic oscillator states in momentum or position 
representation allows a very accurate implementation of overlaps in both representations. Hence, the implementation 
of the Gamow shell-model matrix becomes very similar from one representation to another. 

The Coulomb interaction, not considered in this paper, may however generate difficulties in the momentum repre- 
sentation. Its infinite range character can indeed be treated exactly in the r-representation at the basis level through 
the use of Coulomb wave functions, whereas approximations have to be performed in &;-space calculation as the 
Fourier-Bessel transform of the Coulomb interaction does not exist. The use of the harmonic oscillator expansion 
method may be indeed slowly converging for its low multipoles. This question will have to be answered with Gamow 
shell-model calculations of nuclei close to the proton drip-line. 

VII. CONCLUSION AND FUTURE PERSPECTIVES 

In this paper we have presented a calculational algorithm which can be used to obtain a self-consistent single- 
particle basis in the complex energy plane, starting from a renormalized and realistic nucleon-nucleon interaction. In 
this work we used the simplest possible approximation to the nucleon self-energy, including the Hartree-Fock diagram 
only in order to demonstrate the feasibility of our method. With our approach we studied the single-particle spectra 
of 4 He and 16 0. For 4 He we found that both the P3/2 and the P1/2 states appeared as resonances. Their widths are 
in fair agreement with experiment. For 16 we found the hole states to be largely overbound, while the S1/2 and d§/2 
particle states agree well with the experimental values. The d 3 / 2 states comes out as a resonance, in agreement with 
experiment although our width is larger than the experimental value. Higher order corrections such as 2p — lh and 
lh — 2p contributions may improve the agreement with experiment, in particular the spin-orbit splittings for the hole 
states, and will be included in future self-energy calculations. 

With the Gamow Hartree-Fock single-particle basis, derived from realistic interactions, the problem of representing 
the nucleon-nucleon interaction in the derived basis comes to the fore. As the nucleon-nucleon interaction is typically 
given in momentum space, a transformation from the relative and center of mass frame to the laboratory frame 
involves mathematical functions which are not easy to continue analytically in the complex fc-plane. To that end, we 
investigated whether a method based on an expansion of the nucleon-nucleon interaction as function of a finite set 
of harmonic oscillator functions could be a promising route. This expansion allows for a straightforward calculation 
of matrix elements in the laboratory frame for any Gamow basis. The harmonic oscillator functions are indeed very 
flexible in both position and momentum space, and the analytic continuation of the two-body interaction in the 
complex plane turns out to be very easy to implement from both a theoretical and a numerical point of view. With 
this method, we have shown that for the example of the 6 He and 18 O nuclei, all states converge with a low number of 
harmonic oscillator functions in the expansion of the interaction. This method offers also a practical way of calculating 
higher order diagrams in many-body perturbation theory. There are diagrams which enter for example the definition 
of the self-energy and two or three-body effective interactions. In particular, this method provides a solution to the 
problem of non-orthogonal intermediate particle states, a problem which arises when one uses the vector brackets. 
Utilizing the harmonic oscillator expansion of the interaction, we ensure that all intermediate states are orthogonal 
in all diagrams beyond first order in many-body perturbation theory. However, the renormalization of the Coulomb 
interaction and the two-body center of mass contribution need further considerations. 

The procedure outlined in this work allows for several interesting applications related to the study of weakly bound 
and unbound nuclei along the driplines. Of particular interest is the possibility to apply our approach within the 
framework of the coupled-cluster method, see for example Refs. 0, HE EE US ES H3 • The coupled-cluster approach 
is a promising candidate for the development of practical methods for fully microscopic ab initio studies of nuclei. 
The coupled-cluster methods are capable of providing a precise description of many-particle correlation effects at 
relatively low computer costs, when compared to shell-model or configuration interaction techniques aimed at similar 
accuracies. In approaches such as the coupled cluster, an extension to the complex energy plane is in principle 
possible. Coupled-cluster calculations starting with a HF basis in the complex plane, along with the inclusion of 



16 



realistic interactions, are planned in near future. Then, it might be possible to perform coupled-cluster calculations 
of states with a multi-particle resonant structure starting with a realistic interaction. 
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